Relaxation of Fermionic Excitations in a Strongly Attractive Fermi Gas in an Optical Lattice 
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We theoretically study the relaxation of high energy single particle excitations into molecules in a system of 
attractive fermions in an optical lattice, both in the superfluid and the normal phase. In a system characterized 
by an interaction scale U and a tunneling rate t, we show that the relaxation rate scales as ^ Ct exp(— at/^/t^) 
in the large U/t limit. We obtain explicit expressions for the exponent a, both in the low temperature superfluid 
phase and the high temperature phase with pairing but no coherence between the molecules. We find that 
the relaxation rate decreases both with temperature and deviation of the fermion density from half-filling. We 
show that quasiparticle and phase degrees of freedom are effectively decoupled within experimental timescales 
allowing for observation of ordered states even at high total energy of the system. 

PACS numbers: 



Ultracold atoms on optical lattices 1 1 1 can be used for simu- 
lating strongly interacting quantum many body systems ||2j]4l 
with tunable Hamiltonian parameters like interaction strength. 
Although the main focus of cold atom experiments has been to 
obtain the equilibrium phase diagram of various models, ultra- 
cold atomic systems also provide an unique platform to study 
the intrinsic non-equilibrium dynamics of strongly interacting 
many body systems. They are almost completely decoupled 
from external environment and their low energy- scales lead to 
long non-equilibrium time-scales over which it is possible to 
follow the system without ultra-fast probes. 

The issues of non-equilibrium relaxation dynamics are be- 
coming an important consideration for the state of the art cold 
atom experiments [ 5 -111 as well as recent pump-probe exper- 
iments with electron systems |[T2l . The tunability of Hamilto- 
nian parameters to access strongly interacting regimes is one 
of the central attractive features of cold atoms. However, an 
implied assumption in connecting the results obtained on op- 
tical lattices to the physics of condensed matter systems is that 
the atoms on the optical lattice have achieved thermal equilib- 
rium at low temperatures after tuning the parameters. Hence, 
it is important to understand the relaxation dynamics and as- 
sociated equilibration timescales 1 13 - 20] of these systems. 

The attractive (or negative U) Hubbard model on opti- 
cal lattices is a lattice implementation Ell of BCS-BEC 
crossover (U |22l, which is a paradigm for understanding 
strongly interacting superfluids. At weak coupling, this model 
exhibits BCS superfluidity with an exponentially small critical 
temperature. At strong coupling, the physics is governed by 
formation of tightly bound molecules which undergo Bose- 
Einstein condensation at low temperatures. In this Letter, 
we consider the relaxation dynamics of the attractive Hub- 
bard model in a cubic lattice in the strong coupling limit, 
where most fermions are paired to form molecules. We fo- 
cus on the decay of excess unpaired fermions present in the 
system (either due to an external drive like lattice modulation 
or due to sweeping of the Hamiltonian parameters) to form 
molecules. For these high energy excitations, energy conser- 



vation requirements lead to a very slow decay rate that scales 
super-exponentially with the ratio of the interaction strength to 
the bandwidth of the system. Using a particle-hole transform 
to map this problem to that of spin mediated decay of dou- 
ble occupancies in the repulsive Hubbard model, we compute 
the decay rate both in the low temperature superfluid phase 
and in the high temperature paired phase for arbitrary filling 
fractions in the lattice. We find that the decay rate decreases 
both with temperature and with the deviation of the fermion 
density from half filling on either side. We discuss the impli- 
cations of these results for maintaining adiabaticity during a 
sweep of Hamiltonian parameters. 

We consider the one band attractive Hubbard model for 
fermions on a 3D cubic optical lattice 

^ = E 4^.- - ^ E ^^t^u (1) 

{ij) ^ 

where t is the tunneling matrix and U is the local attraction 
between the fermions. In the strong coupling limit (large 
U/t), the fermions are paired to form tightly bound bosonic 
molecules with a large binding energy ~ U, which undergo 
Bose condensation at a temperature ^ J = At'^ /U, controlled 
by the kinetic energy scale of the molecules. The separation of 
the energy scales U ^ t ^ J allows us to consider two differ- 
ent temperature regimes where most of the fermions are paired 
into molecules: (a) the low temperature (T <C J) superfluid 
phase, where the molecules are Bose condensed and (b) the 
high temperature (T t ^ J) phase where the molecules 
do not have phase coherence. We calculate the decay rate of 
unpaired fermions in these two regimes. 

The simplest process, where two unpaired fermions hop on 
top of each other to form a molecule, is forbidden unless the 
binding energy of the molecule (~ U) is carried off by other 
excitations in the system. There are two different modes of 
excitations where this excess energy can be dumped: (a) ki- 
netic energy of other unpaired fermions, (with a scale ~ t) and 
(b) kinetic energy of the molecules, with a scale ~ J. In this 
paper, we assume the later process is dominant, which limits 
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TABLE I: Equivalence of different quantities under the mapping be- 
tween the attractive and the repulsive Hubbard model 



Attractive Model 


Repulsive Model 


Unpaired Fermions 


Doublon hole pairs 


Binding Energy 


Mott gap 


Deviation from half-filling 


Magnetization 


Superfluid order 


Canted antiferromagnetic order 


K. E. of Molecules 


Superexchange energy 



the density of unpaired fermions to be less than J/t ~ t/U. 
Since the kinetic energy has an energy scale of J = 4t^//7, 
n U/ J ~ U'^ /t^ molecular excitations have to be created 
to release the binding energy. The matrix element for this pro- 
cess, within n^^ order perturbation theory, is given by 



M 



t t 
J2J' 



t 

nJ 



t 



Ct exp 



-a-^\n{U/t) 



(2) 

where we have used n\ = for large n and nJ = U to 
write the final form. The decay rate F ~ thus decreases 
super-exponentially with U/t. 

To obtain a physical picture of the decay process, it is 
instructive to use a particle-hole transformation [23 1 which 
maps the attractive (negative U) Hubbard model to a repul- 
sive (positive U) Hubbard model. The attractive model with 
zero magnetization (equal up and down spin densities) at any 
density is equivalent to the repulsive model at half-filling (one 
particle per site) with a finite magnetization proportional to the 
deviation of the fermion density in the attractive model from 
half-filling, i.e. m = (1/2)(1— p), where p is the fermion den- 
sity in the attractive Hubbard model. Under this transforma- 
tion, molecule formation is mapped to the formation of a Mott 
insulator, and the unpaired fermions are equivalent to the high 
energy double occupancy (doublon) hole excitation, with the 
binding energy of the molecules playing the role of the Mott 
gap. At half-filling, the low energy physics of the repulsive 
Hubbard model reduces to an antiferromagnetic Heisenberg 
model which exhibits a canted antiferromagnetic order in its 
ground state in the presence of finite magnetization. This spin 
ordering is equivalent to the emergence of superfluidity in the 
attractive model with the kinetic energy of the molecules play- 
ing the role of spin wave fluctuations. The equivalence of the 
relevant quantities under the mapping is shown in Table. |l] 
Here, we will use the language of the spin model to look at 
quantitative estimates of the decay timescales. 

Decay in the superfluid phase: The superfluid phase of the 
attractive fermions is represented by the canted antiferromag- 
netic phase for the spins in the repulsive model. As a hole hops 
in the background of a Mott insulator with canted spin order- 
ing (shown in Fig.[TJa)) from a site z to a neighboring site j, it 
pushes back the spin on the site j to the site i. This disrupts the 




FIG. 1: Hopping of a hole in a canted antiferromagnetic back- 
ground, (a) A single hop moves back one spin, creating broken bonds 
shown by dashed lines, (b) and (c) Configurations before and after 
multiple hops of a hole. The solid line denotes the trajectory of the 
hole, while the dashed lines in (c) shows the broken bonds. 



spin texture and creates purely ferromagnetic bonds between 
nearest neighbors, each of which gains an energy of ( Jx/2), 
where x = 1 — 4m^ is proportional to the antiferromagnetic 
component of the spin order. Hopping of the hole along a path 
creates ferromagnetic bonds in the directions transverse to this 
path, thus creating a domain wall in the system, as shown in 
Fig. [TJb) and (c). It is to be noted that hopping of doublons 
in this background also leads to a similar process. In a cubic 
lattice each hop creates z — 2 = 4 broken bonds. So, in order 
to accommodate an energy U and decay, the hole (doublon) 
need to traverse a path of length n = 2U/[{z — 2)Jx]. Within 
Fermi's golden rule, the decay rate is given by 



27r peMn)\MfiWn) 



(3) 



where pex is the density of holes (or equivalently of unpaired 
fermion excitations), and Mfi{n) is the matrix element con- 
necting the initial state with a doublon and a hole to the final 
state with a domain wall of length n given by 



Mfi{n) = 



t 



nl 



2t 



{z - 2)Jx 



t(t/UY 



(4) 



Here y{n), the density of states, is proportional to the number 
of self-avoiding paths of length n connecting the doublon and 
the hole and is given by 



u{n) 



2n 



{z - 2)Jx 



S{n,r)G{r) 



(5) 



where 5'(n, r) is the number of self-avoiding paths of length 
n connecting two points at a distance r and G{r) is the di- 
mensionless doublon-hole pair correlation function. We as- 
sume that the doublons and holes are uncorrected, i.e. the 
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FIG. 2: The decay rate of unpaired fermions in the superfluid phase 
plotted as a function of U /t for different fermion densities: (i) half- 
fiUing (soUd black line) (ii) pf — 0.7 (red dashed line) and (iii) 
pf = 0.5 (blue dotted line). 




FIG. 3: The decay rate of unpaired fermions in the high temperature 
normal state plotted as a function of U/t for different fermion densi- 
ties: (i) half-filling (solid black line) (ii) pf — 0.7 (red dashed line) 
and (iii) pf = 0.5 (blue dotted line). 



pair distribution function is independent of r, or G(r) = 1. 
This assumption is valid in the limit of low density of dou- 
blons, precisely the limit we are interested in. Then v(n) = 
{z-2)Jx ^^^)' where S{n), the total number of self-avoiding 
paths of size n, scales as ~ g'^n^ for large n. The constants 
g = 4.68 1 24 1 and /c = 1/6 [1251 for a cubic lattice has previ- 
ously been computed in the context of polymer physics. Using 
these, finally obtain the relaxation rate of unpaired fermions in 
the attractive Hubbard model 

(6) 

Note that, for the attractive case, x = 2p — is a measure of 
the filling factor which vanishes both at p = (empty band) 
and p = 2 (completely filled band) and attains its maximum 
value at half-filling (p = 1). The low temperature decay rate, 
plotted as a function of U/t for different densities in Fig. [2j 
show the expected super-exponential scaling. As we move 
away from half-filling (with increasing magnetization), the 
energy lost in a single hop decreases, and longer domain walls 
are required to absorb the excess energy, leading to a slower 
decay rate. 

Decay in the High Temperature Normal State: We now con- 
sider the decay of the single particle excitations in the high 
temperature phase (T ^ t J), where the fermions are still 
paired into molecules but there is no superfluidity. In terms 
of the spin model, this regime corresponds to a completely 
spin disordered phase with no spatial or dynamic correlations. 
The single site Hilbert space can be occupied by an t spin 
with probability 1 /2 + m or by a down spin with probability 
1/2 — m (we are working at fixed magnetization). 

As in the low temperature phase, the motion of holes (dou- 
blons) pushes the spins along the trajectory by one site. How- 
ever, the energy lost in a given hop now depends on the con- 
figuration of the neighboring spins along the path. Since there 



is no spatial correlation between the spins, the energy lost 
in each hop can be treated as an independent random vari- 
able which takes the values Jr/2 (r = —4, —3, ...4), with 
the probability P(r) = {a/by^^ Y^^^I^ 4^+^ a'b^-\ 
where nCk is the Binomial co-efficient, a = (1/2 — m)^ 
and 6 = (l/2 + m)^. The mean energy lost in each hop is 0, 
while the variance of the distribution is given by J^x/2. The 
total energy lost in / steps is thus a Gaussian random variable 
with zero mean and a variance lJ'^x/2. 

P{E, I) = —1^ eM-EV{J^lx)] (7) 

Since the doublon needs to lose an energy U to decay, one 
must now average over decay processes from paths of length 
/ > n with the probability distribution P{U^l). The square 
of the matrix element for a process involving / hops is given 
by ^t\tyel){ty el)... {tyef) ^t\2tyiJ^x)K where 
Ci is the energy lost after i steps and, to leading order, we 
have replaced ef by its average value iJ'^x/2. To see the 
scaling in t/U, we note that the Gaussian probability dis- 
tribution P{U^l) would be cut-off at a typical lengthscale 
/ ~ 2/7^/(xJ^), and hence the square of the matrix element 
scales SiSt'^{t/UyK Then, summing over all paths with / > n 

F ^ 2^tVe. 1^ dlgH^P{U, I) (^^^ (8) 

The high temperature decay rate, plotted as a function of U/t 
in Fig. [3] is orders of magnitude smaller than the low temper- 
ature decay rate. To understand this, note that the motion of 
the high energy pair can both excite and de-excite low energy 
modes in the background system. At zero temperature, these 
low energy modes are unoccupied and the motion of the high 
energy pair then leads to excitation of these modes. As tem- 
perature increases, the occupation probability of the low en- 
ergy modes increases and the motion of the high energy pair 
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randomly leads to excitation and de-excitation of these modes. 
Thus the high energy pair loses its energy more efficiently at 
lower temperatures and decays faster. 

Adiabaticity and sweep rates: In cold atom experiments the 
strongly interacting regime of model systems is accessed by 
tuning the Hamiltonian parameters at a finite rate, which is 
limited by the lifetime of the atoms in the trap. The tuning pro- 
cess needs to be adiabatic in order to remain in the interesting 
low temperature regime for the system. Since the microscopic 
relaxation processes determine the timescale for equilibration, 
it is expected that the relaxation timescales, along with experi- 
mental sweep rates, would determine the limits of adiabaticity 
in these experiments. We now make these ideas more precise 
by looking at the constraints due to the slow decay of unpaired 
fermions (or equivalently doublon-hole pairs). 

We consider a system of attractive fermions at low enough 
temperatures so that most of the fermions are paired to form 
molecules. In the large U/t limit, the density of unpaired 
fermions in equilibrium p^^ ^ exp(— /7/T), where T is the 
temperature of the system. We assume an adiabatic sweep of 
U/t dXdi constant rate 7 = {U/t) and try to assess the limits 
where adiabaticity fails. At these low temperatures, almost all 
the entropy comes from the kinetic motion of the molecules; 
so for a constant entropy process, we can assume T/ J = A/4 
oxU /T — U^ / (At^) along the sweep , where A is a constant. 
Now, adiabaticity will be maintained in the regime where 

Pex = -{U/T)pea: = 'JJ^Pex < -T{U/t)pe:, (9) 

where 7 is the experimental sweep rate. As the microscopic 
rate T goes down super-exponentially with U/t, this criterion 
would set an upper limit of {U /t)max{l), which is the maxi- 
mum U /t upto which the system remains adiabatic when the 
parameters are swept with a rate 7. Our analysis shows that it 
is extremely difficult to keep the system fully adiabatic in the 
strong coupling limit when either the tunneling or the interac- 
tion are being changed, as the relaxation timescale of unpaired 
atoms (or doublons for the repulsive case) can be anomalously 
long. Experimentally this long timescale should manifest it- 
self as a saturation in the molecular fraction with the satura- 
tion occurring at smaller values of U/t for faster sweep rates. 

At the same time, if relaxation of unpaired fermions is 
very slow (longer then the timescale of experimental measure- 
ments), then they can be considered as infinitely long lived 
and completely decoupled from other degrees of freedom in 
the system like the phase fluctuations of the superfluid order 
parameter. Similarly, in the repulsive Hubbard model, if the 
goal is to observe antiferromagnetism, one may worry that 
a small number of doublons can release enough energy to de- 
stroy magnetic order. If the doublons are very long lived, there 
will be a long timescale over which one can neglect relaxation 
of doublons and analyze the quasi-equilibrium with "unbreak- 
able" doublons. Thus, within experimental timescales, there 
is an effective spin-charge decoupling which makes it easier 
to observe spin ordering even in the presence of high energy 



charge excitations. The idea of realizing metastable states 
with long lived doublons has also been discussed in the con- 
text of the 7^-paired state in the repulsive Hubbard model 1261 . 

We have studied the decay of unpaired fermions in an at- 
tractive Hubbard model. We have shown that the decay rate 
scales as ~ exp[— a/7^/t^] for large U/t and computed the 
exponent a both at low temperatures (superfluid phase) and 
high temperatures (normal state of molecules). We find the 
decay rate decreases with increase in both temperature and the 
deviation of the fermion density from half-filling. We also dis- 
cussed implications of our analysis for realizing many-body 
states in optical lattices. The downside of the long relaxation 
times is that it is difficult to change parameters of the system 
fully adiabatically. The upside of slow relaxation is that there 
is effective decoupling of different degrees of freedom. So for 
example, one may be able to achieve equilibration of phase 
(magnetic) degrees of freedom, even when there is a finite 
density of unrelaxed single fermions (doublons). 
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